arXiv:math/0112195vl [math.AG] 19 Dec 2001 


Real Quartic Surfaces in MP3 containing 16 Skew 
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Abstract 

In [BN] the authors construct a special complex of degree 20 over M, which for 
an open three dimensional set parametrizes smooth complex surfaces of degree four 
invariant which are Heisenberg invariant, and each member of the family contains 
32 lines but only 16 skew lines. The coordinates of the lines however need not 
be real. For a point I in a Zariski open set of M an algorithm is presented which 
evaluates the real coefficients of I in terms of the K- coordinates of 1. The author 
uses a code in Maple which allows him to construct very explicit examples of real 
smooth Heisenberg invariant Kummer surfaces containing the special configuration 
of 16 real skew lines. An example is presented at the end of the paper. 

MSC (2000): 14Nxx,14Lxx,54Dxx,53Cxx,32Mxx. 


1 Introduction 


In the variables Xo,Xi,X 2 ,Xs,X 4 , x^ consider the equation : 

xl + xl + xl + xl + xl + xl =0 
l/Xg + l/Xi + \lx\ + l/Xg + \lx\ + l/Xg = 0. 


( 1 ) 

( 2 ) 


This is a rational form of the equation which must be written in regular form ! 

Let us consider a solution to these set of equations. For this let Ui = xl and hxing 
M2,M3,M4,M5 then in |[Nie|| we have found an equation of the form Mq + Xuq + X/y = 0 
where X,y are constants depending on the above hxed variables. The solutions can 
however be complex. If we let l = \/—1 and introducing 


Poi =Xo- LXi, P 02 = X2- tX3, Po 3 = X4 - LX3, 

P23 =-{xo + LXi), Pi3 =X2 + LX3, Pi2 =-{x 4 + LX5) 


1 
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if all {pi,j} are real, the pij coordinates will determine a real line I in P 3 contained in 
a real quartic surface, which is invariant under the well known group the Heisenberg 
group of level 2. Such surfaces form a three dimensional parameter space within the thirty 
four dimensional space of quartic surfaces in MP 3 and are therefore very special . Here 
we present two algorithms for writing the equation for the surface and the lines lying 
on them together with one example illustrating the second algorithm. These examples 
are surfaces of degree four embedded as Kummer surfaces associated to abelian surfaces 
of polarization type (1,3) (c.f. |[BN|| and |[Nie|| ). We give an elementary solution to the 
problem stated in the beginning of section ^ which solution is the hrst algorithm. In 
section ^ we introduce the second algorithm . In section ^ we briefly recall the basic 
concepts of Line Geometry and Quartic Surfaces in MP 3 to solve our problem ,use the 
second algorithm in section ^ to construct one example and briefly sketch how the code 
was written in Maple. Additionally we produce the equation of the quartic surface and 
the set of skew lines lying on it. 


2 An elementary Problem 

We consider the following elementary problem which hrst solution was communicated 
by the author by W. Barth in an unpublished manuscript: Find solutions to the set of 
eguations || and ^ such that Xq,X 2 , x^ are real and xi,x^, X 5 are purely imaginary. 

For later use we let 


xi = iyi,X‘i = Ly3,X5 = iy5 

where i is a/—1. By substitution the new form of equation |l| is 

2 I 2 I 2 2 I 2 I 2 

Xq + X2 + X4 = 2/1 + 2/3 + 1/5 

and similarly for the second equation 

1 / X 2 + 1/^2 ^/^2 ^ yy 2 ^ ^/^2 ^ ^/^ 2 _ 

We will normalize equation | as 

2,2,2 1 2,2,2 

Xq + X2 + X4 = 1 = 1/4 + ?/3 + 2/5 


and set 

q = I/xq + I/X 4 + I/X 2 = 1 / 2/1 + 1 / 1/3 + 1 /l/i 
If we introduce the following change of variables: 

_ j xj for i even, 

\ Vi otherwise. 


(3) 

(4) 


we obtain the following 
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Proposition 1 Choose 7^ real such that 


0 < g 4 , gs < 1 

and 

q > ma2;(l/g4 + 4/(l-g4),l/g5 + 4/(l-g5)) 
then there exist xq, X2i X4,, yi, y^, y^ real satisfying equation 
Proof - Considering 

1/go + l/g 2 = g - l/g4 

and setting s = 1 — g4 = go + g2 we obtain the formal eqnation 

P = gog2 = (go + g2)/(g-l/g4), 

= (l 4 s/{qq 4 - 1), 

= g 4 (l - q 4 )/{qq 4 - !)• 

Assnming that s,p are given then the roots of 

X^-sX+p = 0 ( 5 ) 

determine go,g 2 whenever g is given. In fact, go = (s + \/s"^ — 4p)/2, 
g 2 = (s — \/s'^ — 4p) /2 . Similarly if s' = 1 — gs = gi + go we obtain the qnadratic 


eqnation: 


-PX + p' = 0 

and the condition q^ ^ q^ implies that s 7^ s'. Eqnation | has the elementary solntions 

(s ± y/s^ — 4 p)/ 2 . 

A necessary condition that eqnation has different solntions is that 

> Ip = 4sg4/(gg4 — 1). 

By assnmption s > 0 thns s > Iq^/i^qq^ — 1 ) which gives 

g > l/g4 + 4 /s. 

It also follows from the last ineqnality and eqnation ^ that 

P > 0 go,g2 > 0. 

Similarly ifg > l/g5 + 4 /(l —gs) and 0 < gs < 1 then yi^y^axe roots of positive nnmbers. 
o 

Remark 2 For this condition the normalization of equation ^ is not needed. Moreover, 
the condition chosen for convenience q^ = gs is indeed necessary since it implies that 

2,2 2,2 2,2 n 

Xq + Xi = Xr + Xg = 0:2 + X3 = 0 


and it follows from / jAzefl/ , §5 and %6] that the surface is elliptic ruled or singular along 
four nodes. 
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3 The Algorithm 

Once again we fix the normalization of equation ^ in the coordinates g*, which is 1 = 
+ ? 4 - We are to hnd real positive solutions to the non-linear system of polynomial 
equations given in rational form as : 

A = 1/^2 + 1/go + 1/ 

1 = Qo + Q2 + Qi ( 6 ) 

For the following we let uj = (Agg — go(A — 3) -f 1)/A and N = (Ago — 1)(1 — go), 
= iV(iV —4go). Note that these are polynomial expressions in go. We £x in the sequel 
the numerical value for go. 

Proposition 3 The roots of equation ^ for g 2 are: g 2 = (A ± A)/2(Ago — 1) with 
A^ = X'^u}{l/X - go)(l - go). 

Proof - Claim: 

g 2 = (Ago - 1)(1 - go) ± V(-^^o - 1)^(1 - go)^ - 4go(l - go)(Ago - l)/2(Ago - 1) 

To prove the claim from the above given equation one obtains 

go<?4 + g2g4 + go<?2 - Agog2g4 = 1 - go - <?2 - g 4 = 0 
By substituting the value of in the last equation one obtains the quadratic equation 
ql{Xqo - 1) - g 2 (Ago - 1)(1 - go) + go(l - go) = 0 
which solution is as claimed. 

To hnish the proof we write A^ = NM where M = (Ago — 1)(1 — go) — 4go. The last 
equality is equal to —(AgQ -|- go(A — 3) -f-1). o 

Remark 4 The condition 1 = go + g 2 + g 4 still has to be fulfilled which does not follow 
from Proposition |^. 

Remark 5 Let go > 1/A then 

go + g 2 < 1 A — eA > 0 with e G {±1}. 

Proof of the remark An easy computation shows that 

go + g 2 = ((Ago - l)(go + 1) + eA)/2(Ago - 1). 


Hence obtain (Ago — 1)(2 — (go -|- 1)) -|- eA > 0. 

We introduce more notation. Let a = (A — 3)/2A,p = (A — 9)(A — 1)/4A^. 


o 
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Proposition 6 > 0 oo > 0,qo < 1/X or a;<0,l/A<go- 

Proof - This is a consequence of Proposition |^. o 

Lemma 7 a; > 0 go > <^ + a/P or qq < a — . 

Proof - By writing: 

- 9 o(A - 3)/A+1/A, 

= (g2 - 2(A - 3)/2Ago + (A - 3)V4A2) + 1/A - (A - 3)V4A2, 

= {Qo — (A — 3)/2\y — p. 

The cases claimed in the lemma follow inmediately from it. o 

The proof of the following lemmma is analogous to the previous one: 

Lemma 8 a; < 0 go G (a — a + ^/p) . 

o 

Lemma 9 If N > 0 then A — A > 0. If additionally go G (a — a + ^/p) then 
go + g 2 < 1. 

Proof - A^ = —NXuj. We have to evaluate N — \/—NXu. But 

(A - \/-AAa;)(A + V-NXuj) = A(A + Ao^) 

From A + Ao; = 4go > 0 thus A — A > 0. Since A > 0 if and only if go > 1/A the last 
remark follows from Remark |^. o 

We shall need the following easy arithmetic: 

Lemma 10 Assume A > 9 then a > y/p and 1/A < cr — y/p. 

Proof - To prove the hrst inequality: 

A > 0 ^ (A - 3)2 > (A - 9)(A - 1) ^ (A - 3) > y/(A - 9)(A - 1). 

To prove the second inequality: 

2 < (A-3)-y/(A-9)(A-l) 

^ (A-5)>y/(A-9)(A-l) 

(A-5)2 > A2-10A + 9 

and the last inequality is always true. o 

Corollary 11 Let go < 1/A. go + g 2 < 1 —A + A>0 or—A — A>0. 
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Proof - By using the proof of Remark I go + g 2 < 1 -(Ago - l)(go + 1 ) + eA < 
2(1 — Ago) +eA < —N 0 < —N ± A with e G {±1}- ^ 

Proposition 12 Assume A > 9 and Ago < 1 then go + g 2 < 1- 

Proof - Clearly Ago < 1 Ago — 1 < 0 thus A < 0. By lemma to prove a; > 0 

we verify one of the inequalities : go < cr — ^/p which is true by lemma ^ hence 
1/A < a — yp. By Prop. | one of the inequalities is satished thus A^ > 0 and — A +A = 
—N + \/—N\uj > 0 which is one of the inequalities to be satished in Corollary 0 o 

Proposition 13 Assume A > 9 and go G (a — a + y/p) then go > 1/A, go + g 2 < 1- 

Proof - If A > 9 then 1/A < a — ^p by lemma but then go > 1/A. In this case A > 0 
hence by lemma ^ the inequality go + g 2 < 1 is true. o 

Remark 14 The cases A = 1, 9 m the above proposition give no criteria to find go as a 
solution to equation This can be explained as follows: fix the affine plane H = {g = 
go + + g 4 — 1 = 0 } in the defined by the coordinates go, g 2 , Qa and define for each 

A G M>o-' 

/a = ^2^4 + qoQa + qoq2 - Agog2g4 

which is a surface in the A 3 defined by the coordinates go, g 2 , qA- Let C\ = H P {fx = 0}. 

Proposition 15 The linear system of curves {C\ = {fx = g = 0}}agr>o always 
smooth except for A = 1, 9 with singularities: 

SingiCg) = {(1/3,1/3,1/3)}, Sing{Ci) = [Q = (-1,1,1)}. 

Proof - To simplify the computations let dif = df /dqi. Recall that fx = q 2 qA + go ?2 + 
gog 4 -Agog 2 g 4 and-g = l-go-g 2 -g 4 - Note that 9*/ = qk-\qjqk ^01 i,j,k G {0,2,4} 
with i ^ j ^ k and dig = — 1. The Jacobian matrix is 

d2fx dofx difx ^ 

It is of rank less than one if and only if: 

o = (ai-a,)^ = o^,JG{o,2,4}. 

It follows that : (g* — gj)(l — Ag^) = 0 for some i,j, k G {0, 2,4}. Due to the symmetries 
of the indices it is enough to consider the following cases: 

I. 1 - Ag 4 = go - g 4 = g 2 - g 4 = 0 . g 2 = g4 = go- <?4 = 1/A. From 0 = g = 1 - 3go. 
Therefore go = 1/3. Substituting 0 = /a = 3 — A/3. Thus A = 9. 
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II. 1 — Xq 4 = 1 — Xq 2 = q 2 — q 4 = 0. Therefore q 2 = q 4 = I/X. g = l — 2/X — qo = 0. 
Substituting qo = 1— 2/X in f\ = 0 one obtains 0 = f\ = (1 + A — 2)/A = 0. Thus A = 1. 

III. l — Xq 4 = 1 — Xq 2 = 1 — Ago = 0. Therefore go = g 2 = <?4 = 1/A. Hence fx = 2/X^ = 0 
which is impossible. 

IV. 0 = g 2 — go = go — g 4 = g 2 — Q 4 - If follows that go = g 4 = g 2 - Hence g = 1 — 3go = 0. 

Substituting go = 1/3 in /a one obtains /a = 3/9 — A/27 = 0 therefore A = 9. o 

Fix once again the dehned by the coordinates x, y, z. 

Lemma 16 The equation of {C\ = {0 = fx{x,y,z) = 1 — {x + y + 2:)}}AeR>o can be 
written for X = 9 as : 

fg = 9{x‘^y + xy"^) + x + y — 19xy — 

/i = x‘^y + xy"^ + X + y — 2xy — + y"^). 

The proof of the lemma is left as an easy exercise. For the next lemma recall that 
P = (1/3,1/3,1/3) is a singular point for Cg and Q = Sing(Ci). 

Lemma 17 The equation for the tangent cone of Cg (resp. of Ci) passing through 
P(resp. Q of Cl) is (x — 1/3)^ + (x — 1/3)(|/ — 1/3) + (y — 1/3)^ (resp. 4{y — l){x + y)). 
In particular, P (resp. Q) is a non-ordinary double point of Cg at P (resp. Q of Ci). 

Proof - The second partial derivatives at P are given as d^fg = —2 + 18?/, d^yfg = —10 + 
18a; + 18?/, d^fg = -2 + 18a;. Summarizing : dlfg{P) = 4:,dlfg{P) = I,dlyfg{P) = 2. 
The equation for the tangent cone at P is then 

4(a; - 1/3)2 + 4(a; - l/3)(?/ - 1/3) + 4(?/ - 1/3)^. 

The calculation for Ci can be done analogously. o 


Remark 18 A direct computation shows that our fg = 9is irreducible over R . Under 
the linear change of coordinates u = x — 1/3,?; = y — 1/3 the equation for fg = 9 is 
transformed to fg = 9uv{u-\-v)-\-2uv-\-2{u'^-\-v^). Under this linear change of coordinates 
the cubic curve Cg is transformed to a real cubic with isolated singularity at the origin 
which is to be expected from the classification of irreducible cubic curves over the real 
number field. 

4 The 32 lines on the quartic surface. 

We start our dicussion with some well known facts on Line Geometry and Group the¬ 
ory. Fix the three dimensional real projective space RP 3 with coordinates Zg, Zi, Z 2 , Zg. In¬ 
troduce the quartic surface Nf = {f(zo, zi, Z 2 , Z 3 ) = 0 } given by the homogeneous poly¬ 
nomial / of degree four in the variables zg, zi, Z 2 , zg. A line in RPais generated by a two 
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plane in M 4 represented by a two by four matrix 

/ 1 0 * * \ 

0 1 * * y 

(There are in fact 6 ways of choosing the matrix with this property !) In studying 
problems related to the study of the Geometry of hypersurfaces of the grassmanian of 
lines in P 3 , known as Line Geometry, there exist various choices for the coordinates to 
use. In studying the line complex it was typical in the ninteenth and the early twentieth 
century to introduce special cooordinate systems such as elliptic, pentaspherical coordi¬ 
nates (c.f. [[^e], chap. VIII, §130 and chap. XII, §221]. A coordinatefree approach uses 
only linear algebra to characterize properties of the line complex (c.f. El, chap. 6 ]). 
Here we choose the hrst approach due to the symmetries of our problem ( c.f. |[BM|| for 
the relation of this approach to polarized abelian surfaces). More precisely let 

V ^0 ^2 ^3 / 

be a two by four matrix and introduce the following coordinates 

Pij = ZiZj - Zjz'i i,j e {0,1, 2, 3} 

where i ^ j. These are the Thicker coordinates or P-coordinates for short ; In such 
coordinates the condition that a matrix A dehnes a two-plane is given as: 

P01P23 - P02P13 + P03P12 = 0 . 

In terms of multilinear algebra this is nothing else than giving a two form u such that a; A 
a; = 0 (c.f. E3 chap. 1]). In the coordinates {pij} this is a hypersurface of degree two ( 
the Plucker quadric). Let be the subgroup of SL(4, R) spanned by the transformations 
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which satisfy the relations: 

= Ti = id, aiTi = -Tiai 
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for i = 1 , 2 . One obtains a central exact sequence of groups: 


{±1} ^ ^0 


where G' ~ Z 2 ^. The explicite action of on / for a polynomial as above on the P- 
coordinates is induced by the usual linear action induced on the polynomials of degree 
four, in particular 

a{zQZiZ 2 Z^) = zqZiZ 2 Z^ for all a G H^. 


The action of H^, the (unique up to a constant) Schroedinger representation of degree 
four on induces a representation on as given in ||BN| , p. 178]. Introduce the 

following coordinate transformation in 


Xq = Poi - P23, X2 = P02 + Pl3, X 4 = P03 “ Pl2, 

Xi = \/^(poi +P23),X3 = \/^(P02 - P13),X5 = V^(P03 + P 12 ) • 


Note that 


Xo + X^ + X 2 +X 3 +X 4 + X^ = - 2 (poiP 23 “ P02P13 + P03P12) 


(7) 

which is the equation for the Pliicker quadric, which parametrizes the set of lines in 
MP 3 . The importance of the coordinates {xi} called the Klein coordinates , hereafter 
K-coordinates, is that these are eigenfunctions for the action of the Heisenberg group 
on them. These eigenfunctions are useful for studying the conhguration of lines as 
BN|, section 4]. Equation parametrises the set of lines in MP 3 . By inverting the 


m 


transformation induced by equation ^ we obtain the coordinates for the P-coordinates 
in terms of the K-coordinates as presented in the introduction. Using equation ^ the 
statement of the PROBLEM stated in section ^ is formulated equivalently as: 

The line with P-coordinates {pij} defines a real line if and only if xo,X 2 ,X 4 are real and 
xi,X 3 ,X 5 purely imaginary. 

The P-coordinates can be expressed in terms of the {gi}-coordinates as: 


/ Poi \ / \ 

I P 03 = I , 

V Pl 3 / \ ^ j 

/ P02 \ / \ 

P 23 = -(7%- 

V P12 / V -(v^- / 


( 8 ) 


Let us consider a line I C X with coordinates {pij} such that poi 7 ^ 0 ( which is 
equivalent to ^/qo 7 ^ The two points with coordinates 


( \ = ( ^ 0 X y \ 

\ Pb J \0 I u V J 
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in MP 3 are determined by the P-coordinates as: 


U =P02, V =P03, 

y = -Pi3, X = -pi2. 

The line spanned by Pa,Pb will be by dehnition 

Let us consider quartic surfaces which are invariant under the action of if*. We £x 
once again the quartic surface X = {/ = 0} then af = for all a & such 
that Ao- G {±1}. One sees inmediately ( ||Nie| , Prop. 4.1.1 ii)] ) that the set of real 
quartic forms invariant under ii* is a real vector space of dimension five. By fixing such 
a form it depends on the real values Aq, Ai, A 2 , A 3 , A 4 which are the coefficients of that 
form. If a form of degree four contains a line it imposes five linear conditions on the 
coefficients. More precisely, the restriction of / to a line la,b induces hve linear equations 
in the coordinates of the two points on the line and the coefficients Xi for z = 0,..., 4. It 
follows ( ||B1N| , §4.3]) that each line lab with K-coordinates satisfying equation belongs 
to at least one quartic surface dehned by a quartic form invariant under ii*. Therefore 
the hve linear equations are linearly dependent and projectively four of these equations 
determine the set of coefficients {Aj} of the quartic form which depend on the equivalence 
class of la,b modulo the action of ii*. We will assume in section |]that A 4 = 1. The Maple 
program given in section ^ allows us to determine two other points of the line given two 
of the points of the line lab hence by the previous discussion determine the coefficients 
Xi for z = 0,..., 4 of the quartic surface. If an ii* invariant quartic contains a line I then 
the orbit of I under G' contains 16 skew lines. In the K-coordinates {xi} the involution 

: (xo : Xi : X2 : X3 : Xi : X5) ^ (-I/xq : l/xi : I/X2 : I/X3 : I/X4 : l/x^) 


which is well dehned away from the fourfolds {xi = 0} applied to I gives a line V. Writing 
the P-coordinates associated to this line (for this let q' = —1/^/% — ): 


( /pis I = I 1/^94-1/795 I 

Xil'-P'a/ V 1 / 7 ®+ 1 / 7 *./ 

/ 9 '-p7 \ / 1 / 7 ®-1/7® \ 

9''P23 = l/7»-l/7ft 

Ki'-Pa) \- 1 / 7 ® + 1 / 7 ® / 




(9) 


The previous P-coordinates allows us in the Maple program given in section ^ to deter¬ 
mine very explicitely the parametric equation of the transerval lines once given the line 
lab- Using the K-coordinates cuts exactly 10 lines of the Heisenberg orbit of I (c.f. ||BJN| , 
prop. 4.2 b)]). By reasons of degree, I' C X. Hence the orbit of V under G' is contained 
in X. The two orbits of lines that are in X can be grouped as the “even ” lines as those 
having an even number of minus signs in its K-coordinates and the “odd” lines as having 
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an odd number of minus signs in its K-coordinates ( in fact using these K-coordinates 
we have studied in [BN, prop. 4.2] grouptheoretical properties of the conhguration of 
these lines). 

Classically, the set of odd and even lines form a conhguration of type (16io, 16io) on fhe 
quartic surface forming a double 16 ( as is referred in ||BN| , loc. cit]). By [Q a smooth 
quartic surface can contain at most 48 lines so if it contains at least 32 lines it must 
contain exactly 32 lines. 


5 A brief description of the program 


The program written in Maple IV release 4 dehnes the following local variables: 

R, S', gi, go, iV, S'gi, S'g 2 , g 2 , gs, g 4 , <? 5 , iVi, rr, ss, zz, com, ww, m, n, K, quarn 

and the global variable d. It uses lemma ^ given in section ^ The initial value for the 
program is A ( in the program this value is d ). It then evaluates go and {R, S) using 
prop. 1^, this calculation is performed in the procedures Var, Vas. In order to evaluate 
the variable g 2 one needs to introduce the variables M, N in terms of d and go and hnally 
evaluate S'gi = y/MN. The algorithm then evaluates the positive root of g 2 in terms of 
N, S'gi and dqo — 1. It follows that g 4 = 1 — (go + g 2 )- Analogously, if one gives the value 
for gi the program evaluates gs, gs in exactly the same way. 

In order to draw the lines one evaluates the parametric equation for the lines. The 
procedure to evaluate them in the program has been given in section We introduce the 
variables rr, ss in terms of g 4 , gs, g 2 , ga- The last variables are used to give the parametric 
equation of the line I as given by equation The orbit of I under id* is evaluated using 
the procedure grab The parametric equation of the transversal to I is evaluated using 
the variables m, n. For these one uses the variables rr, zz, ww and uses equation IT. The 


coefficients for the quartic surface are evaluated by means of a matrix array using the 
elementary method of the previous section. The equation of the quartic surface in a hxed 
solid herafter named the “clipping solid” is evaluated in a procedure and saved as the 
variable quarn. The procedure graf used to draw the surface (resp. the lines lying on 
this surface) is based on the command (resp. spacecurve) display3d of the library plots 
of Maple. 


6 An Example 

To illustrate the above theory we will introduce a typical example. The value for this 
example was chosen from a series of 43 test examples for d G [9.3,18.00] on very 
different intervals, on rectangular grids of 30^ = 27, 000 points (using the command 
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grid). The program was run on a Samsung pentium 11. The average compilation time 
per image for these examples averaged to three minutes. We introduce the value for d 
which in this case is 18.00 (this is the upper bound!). The necessary condition for qq 
is given in lemma For practical purposes we introduce two small routines written in 
Maple to evaluate qo,qi in terms of d. For this example, go = 0.4168, gi = 0.1713. Us¬ 
ing the value for go and two other procedures written on the program one obtains 
the values M = 2.124999680, iV = 3.792199680. Using S'gi = 2.83874 one evaluates 
g 2 = 0.5101, g 4 = 0.0731. The values for Mi, Ni, Sq 2 ,q 3 ,q 5 ,rr, ss, zz,m,ww,n are 
(respectively) 1.041313580,1.762513580,1.3408, 0.7364, 0.0923, -0.0334,1.5724, 0.14393, 
—0.10242,0.57418, —0.059153. Using the elementary theory described in section one 
obtains a four by four matrix A with entries polynomials of degree four in the variables 
X, y, u, V and a four by one vector such that the extended matrix M = {A, h) gives the 
equation M ■ A = 0 meaning that the quartic surface with vector coefficients A contains 
a line. Written in non-homogeneous form one obtains : A ■ X = b. The procedure mpoly 
substitutes the variables x, y, u, v for the obtained values rr, ss, zz, ww. For this example 
A = (—0.366, —1.44, 0.614, 0.526). We have chosen the clipping solid as {x+v+y+u = 0}. 
In this solid the equation for the quartic polynomial is given by: 

quarn = 0.320x'^ + 0.496g^ -I—3.612 m^ 

-2.796(a;V + a:V + |/V) + uy{-u'^7.224: + 0.992y^) 

+xy{9.Q‘i9x^ + y^9.992) + xu{x^9.Q‘i9 — .221) 

Txyu{—yb.93Q — 14.152n — a;6.288). 

In the image produced by the procedure graf of the Maple program in all the examples, 
the line I has been drawn in red. The remaining colors for the disjoint lines are: blue, 
yellow, sienna, cyan, khaki, pink, turquoise, aquamarine, magenta, plum, violet, braun, 
green, navy, gold. The ten transversals to I have been drawn in color grey. In the plot3d 
command of Maple one has to specify a pair of values {u, n)for the orientation i.e. the 
direction from which the object is to be viewed. For the same group of test examples 
described in the interval mentioned at the beginning of this paragraph we tested pairs 
u G (0,45) and v G (0, 90) showing in each image at most 7 transversals, with no further 
improvement. We illustrate an example at the end of the references for the surface defined 
by the above given polynomial with a grid of 27,000 points allowing us to see seven of 
the transversal lines. It is for d = 18.00 and at {u,v) = (20,36). 
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Figure 1; The quartic surface for d = 18.00 at (20,36). 



